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Abstract 

New results for the double spin asymmetry A [ and the proton longitudinal spin structure function gj 
are presented. They were obtained by the COMPASS collaboration using polarised 200 GeV muons 
scattered off a longitudinally polarised NH 3 target. The data were collected in 2011 and complement 
those recorded in 2007 at 160 GeV, in particular at lower values of x. They improve the statistical 
precision of gj (x) by about a factor of two in the region x < 0.02. A next-to-leading order QCD fit 
to the gi world data is performed. It leads to a new determination of the quark spin contribution to 
the nucleon spin, AE ranging from 0.26 to 0.36, and to a re-evaluation of the first moment of g' 1 . The 
uncertainty of AE is mostly due to the large uncertainty in the present determinations of the gluon 
helicity distribution. A new evaluation of the Bjorken sum rule based on the COMPASS results for 
the non-singlet structure function g^ s (x. Q 2 ) yields as ratio of the axial and vector coupling constants 
kA /gv = 1.22 ±0.05 (stat.) ±0.10 (syst.), which validates the sum rule to an accuracy of about 9%. 
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1 Introduction 


The determination of the longitudinal spin structure of the nucleon became one of the important issues 
in particle physics after the surprising EMC result that the quark contribution to the nucleon spin is very 
small or even vanishing flj. The present knowledge on the longitudinal spin structure function of the 
proton, gj, originates from measurements of the asymmetry A p in polarised lepton nucleon scattering. 
In all these experiments, longitudinally polarised high-energy leptons were scattered off longitudinally 
polarised nucleon or nuclear targets. At SLAC and JLab electron beams were used, electron and positron 
beams at DESY and muon beams at CERN. Details on the performance of these experiments and a 
collection of their results can be found e.g. in Rcf.[2j 

In this Letter, we report on new results from the COMPASS experiment at CERN. By measuring A p , we 
obtain results on gj in the deep inelastic scattering (DIS) region. They cover the range from 1 (GeV/c ) 2 
to 190 (GeV/c ) 2 in the photon virtuality Q 2 and from 0.0025 to 0.7 in the Bjorken scaling variable x. 
The new data, which were collected in 2011 at a beam energy of 200 GeV, complement earlier data taken 
in 2007 at 160 GeV that covered the range 0.004 < x < 0.7 (3j. In the newly explored low-x region, our 
results significantly improve the statistical precision of gj and thereby allow us to decrease the low-x 
extrapolation uncertainty in the determination of first moments. 

In the following section, the COMPASS experiment is briefly described. The data selection procedure is 
presented in Section 3 and the method of asymmetry calculation in Section 4. The results on A p (x, Q 2 ) 
and gj (x. Q 2 ) are given in Section 5. A new next-to-leading order (NLO) QCD fit to the existing nucleon 
g\ data in the region Q 2 > I (GeV/ c ) 2 is described in Section 6 . Section 7 deals with the determination of 
first moments of gj and the evaluation of the Bjorken sum rule using COMPASS data only. Conclusions 
are given in Section 8 . 


2 Experimental setup 

The measurements were performed with the COMPASS setup at the M2 beam line of the CERN SPS. 
The data presented in this Letter correspond to an integrated luminosity of 0.52 fb A beam of positive 
muons was used with an intensity of 10 7 s " 1 in a 10 s long spill every 40 s. The nominal beam momentum 
was 200 GeV/c with a spread of 5%. The beam was naturally polarised with an average polarisation 
Pb = 0.83, which is known with a precision of 0.04. Momentum and trajectory of each incoming particle 
were measured in a set of scintillator hodoscopes, scintillating fibre and silicon detectors. The beam 
was impinging on a solid-state ammonia (NH 3 ) target that provides longitudinally polarised protons. 
The three protons in ammonia were polarised up to |Pr| ~ 0.9 by dynamic nuclear polarisation with 
microwaves. For this purpose, the target was placed inside a large-aperture superconducting solenoid 
with a field of 2.5 T and cooled to 60mK by a mixture of liquid 3 He and 4 He. The target material was 
contained in three cylindrical cells with a diameter of 4 cm, which had then - axes along the beam line and 
were separated by a distance of 5 cm. The outer cells with a length of 30cm were oppositely polarised 
to the central one, which was 60cm long. In order to compensate for acceptance differences between the 
cells, the polarisation was regularly reversed by rotation of the magnetic field direction. In order to guard 
against unknown systematic effects, once during the data taking period the direction of the polarisation 
relative to the magnetic field was reversed by exchanging the microwave frequencies applied to the cells. 
Ten NMR coils surrounding the target material allowed for a measurement of Pr with a precision of 
0.032 for both signs of the polarisation. The typical dilution due to unpolarisable material in the target 
amounts to about 0.15. 

The experimental setup allowed for the measurement of scattered muons and produced hadrons. These 
particles were detected in a two-stage, open forward spectrometer with large acceptance in momentum 
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and angle. Each spectrometer stage consisted of a dipole magnet surrounded by tracking detectors. 
Scintillating fibre detectors and pixel GEM detectors in the beam region were supplemented with Mi- 
cromegas and GEM detectors close to the beam and MWPCs, drift chambers and straw detectors that 
covered the large outer areas. Scattered muons were identified in sets of drift-tube planes located behind 
iron and concrete absorbers in the first and second stages. Particle identification with the RICH detector 
or calorimeters is not used in this measurement. The ‘inclusive triggers’ were based on a combination 
of hodoscope signals for the scattered muons, while for ‘semi-inclusive’ triggers an energy deposit of 
hadron tracks in one of the calorimeters was required, optionally in coincidence with an inclusive trigger. 
A detailed description of the experimental setup can be found in Ref.j4| 


3 Data selection 

The selected events are required to contain a reconstructed incoming muon, a scattered muon and an 
interaction vertex. The measured incident muon momentum has to be in the range 185 GeV /c < Pb < 
215 GeV jc. In order to equalise the beam flux through all target cells, the extrapolated beam track is 
required to pass all of them. The measured longitudinal position of the vertex allows us to identify 
the target cell in which the scattering occurred. The radial distance of the vertex from the beam axis 
is required to be less than 1.9 cm, by which the contribution of unpolarised material is minimised. All 
physics triggers, inclusive and semi-inclusive ones, are included in this analysis. In order to be attributed 
to the scattered muon, a track is required to pass more than 30 radiation lengths of material and it has 
to point to the hodoscopes that have triggered the event. In order to select the region of deep inelastic 
scattering, only events with photon virtuality Q 2 > 1 (GeV/c ) 2 are selected. In addition, the relative 
muon energy transfer, y, is required to be between 0.1 and 0.9. Here, the lower limit removes events that 
are difficult to reconstruct, while the upper limit removes the region that is dominated by radiative events. 
These kinematic constraints lead to the range 0.0025 < x < 0.7 and to a minimum mass squared of the 
hadronic final state, W 2 , of 12 (GeV/c 2 ) 2 . After all selections, the final sample consists of 77 million 
events. The selected sample is dominated by inclusive triggers that contribute 84% to the total number of 
triggers. The semi-inclusive triggers mainly contribute to the high-x region, where they amount to about 
half of the triggers. In the high-<2 2 region the semi-inclusive triggers dominate. 


4 Asymmetry calculation 


The asymmetry between the cross sections for antiparallel (t|) and parallel (tt) orientations of the lon¬ 
gitudinal spins of incoming muon and target proton is written as 


A 


p _ 
LL — 


gti - qtt 

(jtf + (jh ’ 


( 1 ) 


This asymmetry is related to the longitudinal and transverse spin asymmetries A p { and A''], respectively, 
for virtual-photon absorption by the proton: 


Al L = D(A p l +7 1 Al). 

The factors 

y(l — y — T^y 2 /4 — y 2 nr / Q 2 ) 

71 (1 + 'Y 1 y/2){\-y/2)-y 2 m 2 /Q 2 

and 

n= _ y((l + y 2 y/2)(2 —y)-2y 2 m 2 /g 2 ) _ 

y 2 (l - 2m 2 /<2 2 )(l + y 2 ) + 2(1 +R)( 1 -y- 7 2 y 2 /4) 


( 2 ) 

(3) 

(4) 
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depend on the event kinematics, with y = 2Mxj The virtual-photon depolarisation factor D depends 

also on the ratio R = <7l/ot, where oy (ctt) is the cross section for the absoiption of a longitudinally 
(transversely) polarised virtual photon by a proton. The asymmetry A^ is defined as 


A p = g l/2 ~ a 3/2 
* /2 + °3/2 


(5) 


where 01 / 2 ( 03 / 2 ) is the absorption cross section of a transversely polarised virtual photon by a pro¬ 
ton with total spin projection \ (|) in the photon direction. Since both ?] and A p [5] are small in the 
COMPASS kinematic region, A^ ~ A^/D and the longitudinal spin structure function is given by 


Si = 


F p 

"2 


2x (1 +R) 


A p 


( 6 ) 


where F p denotes the spin-independent structure function of the proton. 


The number of events, Nj, collected from each target cell before and after reversal of the target polarisa¬ 
tion is related to the spin-independent cross section o = C] / 2 + O 3 / 2 and to the asymmetry A p as 

Nj = aityiniCi{\ + PsPrfDA p ) , i = ol,cl,o2,c2 . (7) 

Here, a,- is the acceptance, </), the incoming muon flux, n, the number of target nucleons and / the dilution 
factor, while Pq and P\ were akeady introduced in Section 2. Events from the outer target cell are 
summed, thus the four relations of Eq. (j7]) corresponding to the two sets of target cells (outer, o and 
central, c) and the two spin orientations (1 and 2) result in a second-order equation in A’j’ for the ratio 
(.N tA N c 2 )/(N cl N o2 ). Fluxes and acceptances cancel in this equation, if the ratio of acceptances for the 
two sets of cells is the same before and after the magnetic field rotation [6j. In order to minimise the 
statistical uncertainty, all quantities used in the asymmetry calculation are evaluated event by event with 
the weight factor [6j 

w = PsfD. ( 8 ) 

The polarisation of the incoming muons as a function of the beam momentum is obtained from a 
parametrisation based on a Monte Carlo simulation of the beam line. The effective dilution factor / 
is given by the ratio of the total cross section for muons on polarisable protons to the one on all nuclei 
in the target, whereby their measured composition is taken into account. It is modified by a correction 
factor that accounts for the dilution due to radiative events on unpolarised protons [7], The target polari¬ 
sation is not included in the event weight, because it may change in time and generate false asymmetries. 
The obtained asymmetries are corrected for spin-dependent radiative effects according to Ref. [ 8 ] and for 
the 14 N polarisation as described in Refs. 3,9. It has been checked that the use of semi-inclusive triggers 
does not bias the determination of . 


Systematic uncertainties are calculated taking into account multiplicative and additive contributions to 
A p . Multiplicative contributions originate from the uncertainties of the target polarisation, the beam po¬ 
larisation, the dilution factor (mainly due to the uncertainty of R ) and the depolarisation factor. When 
added in quadrature, these uncertainties result in a total uncertainty AA 1 ™ 11 of 0.07A 1 j ) . They are shown 
in Table [T] which also shows the additive contributions. The largest additive contribution to the sys¬ 
tematic uncertainty is the one from possible false asymmetries. Its size is estimated with two different 
approaches. In the first approach, the central target cell is artificially divided into two consecutive 30cm 
long parts. Combining these two 30cm long targets with the outer ones with the same polarisation, two 
independent false asymmetries are formed. Both are found to be consistent with zero. In order to check 
for time-dependent effects, in the second approach the data sample is divided into sub-samples each 
consisting of periods of stable data taking with both field directions for each target cell. The results for 
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Table 1: Contributions to the systematic uncertainty on A’j’ with multiplicative (top) and additive (bottom) compo¬ 
nents. 


Beam polarisation 

A Pb/Pb 

5% 

Target polarisation 

APf / Py 

3.5% 

Depolarisation factor 

A D{R)/D(R) 

2.0 - 3.0% 

Dilution factor 

A/// 

2% 

Total 

AA“ ult 

~0.07Ai 

False asymmetry 

Afalse 

< 0.84 ■ a stat 

Transverse asymmetry 

v • a 2 

< 10~ 2 

Radiative corrections 

Af c 

O 

1 

4^ 

1 

O 

1 

U> 


Aj obtained from these sub-samples are compared by using the method of “pulls” [10]. No significant 
broadening of pull distributions is observed. These pulls are used to set an upper limit on the system¬ 
atic uncertainty due to false asymmetries Aj alse . Depending on the x-bin, values between 0.4a stat and 
0.84-a stat are obtained. Further additive corrections originate from neglecting A 2 and from the uncer¬ 
tainty in the correction A^ c to the asymmetry A \, which is due to spin-dependent radiative effects. 


5 Results on A 1 , and g 1 . 


The data arc analysed in terms of A\ and g\ as a function of x and Q 2 . The x dependence of A 1 ^ aver¬ 
aged over Q 2 in each x bin is shown in Fig. [I] together with the previous COMPASS results obtained at 
160GeV |3j and with results from other experiments ||T] [TTj - jT4| including those by SMC at 190GeV |T5 1. 
The bands at the bottom represent the systematic uncertainties of the COMPASS results as discussed in 
Section [4] The new data improve the statistical precision at least by a factor of two in the low-x region, 
which is covered by the SMC and COMPASS measurements only. The good agreement between all ex¬ 
perimental results reflects the weak Q 2 dependence of A p . This is also illustrated in Fig.[2j which shows 
A 1 , 1 as a function of Q 2 in sixteen intervals of x for the COMPASS data sets at 160GeV and 200GeV. In 
none of the x bins, a significant Q 2 dependence is observed. The numerical values of A p (x,<2 2 ) obtained 
at 200 GeV are given in the Appendix. 


The longitudinal spin structure function gj is calculated from A p using Eq. (ml), the F 2 P parametrisation 
from Ref. [15] and the ratio R from Ref. 16 The new results are shown in FigTj3] at the measured values 
of Q 2 in comparison with the previous COMPASS results obtained at 160GeV and with SMC results 
at 190 GeV. Compared to the SMC experiment, the present systematic uncertainties are larger due to a 
more realistic estimate of false asymmetries, which is based on real events. 


The world data on g p as a function of Q 2 for various x are shown in Fig. [ 4 ] The data cover about two 
decades in x and in Q 2 for most of the x range, except for x < 0.02, where the Q 2 range is much more 
limited. The new data improve the kinematic coverage in the region of high Q 2 and low x values, which 
gives a better lever arm for the determination of quark and gluon polarisations from the DGLAP evolution 
equations. In addition, the extension of measurements to lower values of x is important to better constrain 
the value of the first moment of g p . 


6 NLO QCD fit of gi world data 

We performed a new NLO QCD fit of the spin-dependent structure function g 1 in the DIS region, 
Q 2 > 1 (GeV/c) 2 , considering all available proton, deuteron and 3 He data. The fit is performed in 
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Figure 1: The asymmetry Aj as a function of x at the measured values of Q 2 as obtained from the COMPASS data 
at 200GeV. The new data are compared to the COMPASS results obtained at 160GeV 0 and to the other world 
data (EMC (T), CLAS (TT), HERMES (T2), E143 Q3), E155 Q3), SMC Q3J). The bands at the bottom indicate 
the systematic uncertainties of the COMPASS data at 160GeV (upper band) and 200GeV (lower band). 
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the MS renormalisation and factorisation scheme. For the fit, the same program is used as in Ref. 

The region W 2 < 10 (GeV/c 2 ) 2 is excluded as it was in 
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which was derived from program 2 in Ref. 
recent analyses [ 18]. Note that the impact of higher-twist effects when using a smaller W 1 cut is con¬ 
sidered in Ref. 19 The total number of data points used in the tit is 495 (see Table |2]), the number of 
COMPASS data points is 138. 


The neutron structure function gf is extracted from the ’He data, while the nucleon structure function 
is obtained as 


£(x,Q i ) = 


1 


1 - 1.5 Wd 


gf(x,Q 2 


(9) 


where cud is a correction for the D-wave state in the deuteron, cod = 0.05 ± 0.01 |24|, and the deuteron 
structure function gj is given per nucleon. The quark singlet distribution Aq s (x), the quark non-singlet 
distributions Aqj(x) and Aq&(x), as well as the gluon helicity distribution A g(x), which appear in the 
NLO expressions for gj, gf and g 1 ^ (see e.g. Ref. 15), arc parametrised at a reference scale Q q as follows: 


a f k {x) = m 


x ak (1 — x )& (1 + YkX ) 
fd x Uk (1 — x)& (1 + YkX )dx 


( 10 ) 


Here, A fk(x) (k = 1...4) represents Aq s (x), Aqj,{x), Aq&(x) and Ag(x) and rj/ ; is the first moment of 
A fk (x) at the reference scale. The moments of Ar /3 and A qg are fixed at any scale by the baryon decay 
constants (F+D) and (3F—D), respectively, assuming SU(2)f and SU(3)f flavour symmetries. The impact 
of releasing these conditions is investigated and included in the systematic uncertainty. The coefficients 
Yk are fixed to zero for the two non-singlet distributions. The exponent /3o, which is not well determined 
from the data, is fixed to 3.0225 ['25} and the uncertainty from the introduced bias is included in the final 
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Figure 2: The asymmetry A, as a function of Q 2 in bins of x obtained from the 200GeV (red squares) and 160GeV 
(blue circles) COMPASS data. The band at the bottom indicates the systematic uncertainty for the 200GeV data. 


0.06 


0.04 


0.02 


0 


10' 2 10' 1 1 

x 

Figure 3: The spin-dependent structure function rgj at the measured values of Q 2 as a function of x. The COM¬ 
PASS data at 200GeV (red squares) are compared to the results at 160GeV (blue circles) and to the SMC results at 
190GeV (green crosses) for Q 2 > 1 (GeV/c) 2 . The bands from top to bottom indicate the systematic uncertainties 
for SMC 190GeV, COMPASS 200GeV and COMPASS 160GeV. 
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Figure 4: World data on the spin-dependent structure function as a function of Q 2 for various values of x with 
all COMPASS data in red (full circles: 160GeV, full squares: 200GeV). The lines represent the Q 2 dependence 
for each value of x, as determined from a NLO QCD fit (see Section [6]). The dashed ranges represent the region 
with W 2 < 10 (GeV/c 2 ) 2 . 


uncertainty. This leaves 11 free parameters in the fitted parton distributions. The expression for y 2 of 
the fit consists of three terms, 


N a 


* 2 = L 


n—l 




(=1 


<yV„o, 


+ 




+ Xi 


positivity * 


(ii) 


Only statistical uncertainties of the data are taken into account in a,. The normalisation factors 
of each data set n are allowed to vary taking into account the normalisation uncertainties 5./F. If the 
latter arc unavailable, they are estimated as quadratic sums of the uncertainties of the beam and target 
polarisations. The fitted normalisations are found to be consistent with unity, except for the E155 proton 


data where the normalisation is higher, albeit compatible with the value quoted in Ref. 14 


In order to keep the parameters within their physical ranges, the polarised PDFs are calculated at ev¬ 
ery iteration of the fit and required to satisfy the positivity conditions |A q(x) + Aq(x)\ < q(x) + q(x') 
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and |Ag(x)| < g(x) at Q 2 = I (GeV/c) 2 , which is accomplished by the £ 2 ositivity term in Eq. (fill). This 
procedure leads to asymmetric values of the parameter uncertainties when the fitted value is close to the 
allowed limit. The unpolarised PDFs and the corresponding value of the strong coupling constant oc s (Q 2 ) 


are taken from the MSTW parametrisation [251. The impact of the choice of PDFs is evaluated by using 


the MRST distributions [261 for comparison. 


In order to investigate the sensitivity of the parametrisation of the polarised PDFs to the functional forms, 
the fit is performed for several sets of functional shapes. These shapes do or do not include the ys and 
y g parameters of Eq. (10) and are defined at reference scales ranging from 1 (GeV /c) 2 to 63 (GeV /c) 2 . It 
is observed |27| that mainly two sets of functional shapes are needed to span almost entirely the range 
of the possible A q s (x) and A g(x) distributions allowed by the data. These two sets of functional forms 
yield two extreme solutions for Ag(x). For y g = ys = 0 (y g = 0 and ys / 0) a negative (positive) solution 
for A g(x) is obtained. Both solutions are parametrised at Qq = I (GeV/ c) 2 and lead to similar values of 
the reduced % 2 of the tits of about 1.05/d.o.f. Changes in the fit result that originate from using other 
(converging) functional forms are included in the systematic uncertainty. 


The obtained distributions are presented in Fig. [5] The dark error bands seen in this figure stem from gen¬ 
erating several sets of gi pseudo-data, which are obtained by randomising the measured g i values using 
their statistical uncertainties according to a normal distribution. This corresponds to a one-standard- 
deviation accuracy of the extracted parton distributions. A thorough analysis of systematic uncertainties 
of the fitting procedure is performed. The most important source is the freedom in the choice of the 
functional forms for A q s (x) and Ag(x). Further uncertainties arise from the uncertainty in the value of 
oc s (Q 2 ) and from effects of SU(2)y and SU(3)/ symmetry breaking. The systematic uncertainties are 
represented by the light bands overlaying the dark ones in Fig. [5] For both sets of functional forms dis¬ 
cussed above, As(x) stays negative. It is different from zero for x > 0.001 as are A d(x) and A u(x). The 
singlet distribution A q s (x) is compatible with zero for x < 0.07. 


The inclusion of systematic uncertainties in the fit leads to much larger spreads in the first moments as 
compared to those obtained by only propagating statistical uncertainties (see Table [3]). In this table, A£ 
denotes the first moment of the singlet distribution. Note that the first moments of An + Au, Ad + Ad 
and As + As arc not independent, since the first moments of the non-singlet distributions are fixed by the 
decay constants F and D at every value of Q 2 . The large uncertainty in Ag(x), which is mainly due to the 
freedom in the choice of its functional form, does however not allow to determine the first moment of 
A g(x) from the available inclusive data only. 

The fitted gj and g) 1 distributions at Q 1 = 3 (GeV/c) 2 are shown in Fig. [6] together with the data evolved 
to the same scale. The two curves correspond to the two extreme functional forms discussed above, 
which lead to either a positive or a negative Ag(jt). The dark bands represent the statistical uncertainties 
associated with each curve and the light bands represent the total systematic and statistical uncertainties 
added in quadrature. The values for gj are positive in the whole measured region down to x = 0.0025, 
while gj is consistent with zero at low x. 


7 First moments of g\ from COMPASS data and Bjorken Sum Rule 

The new data on g^ together with the new QCD fit allow a more precise determination of the first mo¬ 
ments ri(2 2 ) = / 0 gi(x,Q 2 )dx of the proton, neutron and non-singlet spin structure functions using 
COMPASS data only. The latter one is defined as 

gf (x,Q 2 )=g p l (x,Q 2 )-g D 1 (x,Q 2 ) = 2[g p l (x,Q 2 )-g^(x,Q 2 )]. (12) 
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Figure 5: Results of the QCD fits to g\ world data at Q 2 = 3 (GeV/c) 2 for the two sets of functional shapes as 
discussed in the text. Top: singlet ;rA< 7 S (;r) and gluon distribution vAg(x). Bottom: distributions of x [Ag(x) +A^(jc)] 
for different flavours (m, d and ,v). Continuous lines correspond to the fit with ys = 0, long dashed lines to the one 
with ys / 0. The dark bands represent the statistical uncertainties, only. The light bands, which overlay the dark 
ones, represent the systematic uncertainties. 




Figure 6: Results of the QCD fits to gj (left) and gf (right) world data at Q 2 - 3 (GeV/c) 2 as functions of x. 
The curves correspond to the two sets of functional shapes as discussed in the text. The dark bands represent the 
statistical uncertainties associated with each curve and the light bands, which overlay the dark ones, represent the 
systematic uncertainties. 
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Table 3: Value ranges of first moments of quark distributions, as obtained from the QCD fit when taking into 
account both statistical and systematic uncertainties, as detailed in the text. 


First moment 

Value range at Q 2 = 3 (GeV/c) 2 

AT 

Au ± Au 

Ad ± Ad 
As± As 

[ 0.26, 0.36 ] 

[ 0.82, 0.85 ] 

[ -0.45 , -0.42 ] 

[ -0.11 , -0.08 ] 


Table 4: Contribution to the first moments of g i at Q 1 = 3 (GeV/c) 2 with statistical uncertainties from the COM¬ 
PASS data. Limits in parentheses are applied for the calculation of Y ). The uncertainties of the extrapolations are 
negligible. 


x range 


■pN 

1 1 

0 - 0.0025 (0.004) 
0.0025 (0.004) - 0.7 

0.7 - 1.0 

0.002 

0.134 ±0.003 

0.003 

0.000 

0.047 ±0.003 
0.001 


Table 5: First moments of gi at Q 2 = 3 (GeV/c) 2 using COMPASS data only. 



r 1 

srf at 

5Tf l 

Sry o1 

Proton 

0.139 

±0.003 

±0.009 

±0.005 

Nucleon 

0.049 

±0.003 

±0.004 

±0.004 

Neutron 

-0.041 

±0.006 

±0.011 

±0.005 


The integral r^ s (<2 2 ) at a given value of Q 2 is connected to the ratio g\/g\ of the axial and vector 
coupling constants via the fundamental Bjorken sum rule 


rf(e 2 ) = fsY(^)<h=\ 


gA 

gw 


cf(G 2 ), 


(13) 


where Cp s (<2~) is the non-singlet coefficient function that is known [28] up to the third order in a s (Q~) 
in perturbative QCD. 


Due to small differences in the kinematics of the data sets, all points of the three COMPASS gi data sets 
(Table [ 5 ]) are evolved to the Q 2 value of the 160 GeV proton data. A weighted average of the 160 GeV 
and 200 GeV proton data is performed and the points at different values of Q 2 and the same value of * 
are merged. 

For the determination of 0, and Fj 1 . the values of g^ and gj 1 are evolved to Q 2 = 3 (GeV/ c) 2 and the 
integrals are calculated in the measured ranges of x. In order to obtain the full moments, the QCD fit is 
used to evaluate the extrapolation to ^ = 1 and x = 0 (see Table [4]) . The moment Tj is calculated using 
gj = 2g^ — gj. The systematic uncertainties of the moments include the uncertainties of Pq, Pj-, f and 
D. In addition, the uncertainties from the QCD evolution and those from the extrapolation are obtained 
using the uncertainties given in Section [6] The full moments are given in Tabled Note that also is 
updated compared to Ref. fT7| using the new QCD fit. 

For the evaluation of the Bjorken sum rule, the procedure is slightly modified. Before evolving from the 
measured Q 2 to Q 2 = 3 (GeV/c) 2 , g'p is calculated from the proton and deuteron gi data. Since there 
is no measured COMPASS value of gj 1 corresponding to the new gj point at x = 0.0036, the value of 
gj from the NLO QCD fit is used in this case. The fit of g^ is performed with the same program as 
discussed in the previous section but fitting only the non-singlet distribution Aq^ix. Q 2 ). The parameters 
of this tit are given in Table [6] and a comparison of the fitted distribution with the data points is shown in 
Fig.[7J The error band is obtained with the same method as described in the previous section. 
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Table 6: Results of the fit of Aq 3 (x) at Q q 
1(G eV/c) 2 . 


Param. 

Value 

V3 

1.24 ± 0.06 

a 3 

- 0.11 ± 0.08 

ft 

99 + 0.5 
z - z - 0.4 

x 2 /ndf 

7.9/13 


Table 7: First moment T^ s at Q 2 = 3 (GeV /c) 2 from the 
COMPASS data with statistical uncertainties. Contribu¬ 
tions from the unmeasured regions are estimated from the 
NLO fit to £|'" s . The uncertainty is determined using the 
error band shown in Fig. [7] 


x range 

pNS 

1 1 

0 - 0.0025 
0.0025 - 0.7 

0.7 - 1.0 

0.006 ± 0.001 
0.170±0.008 

0.005 ±0.002 

0 - 1 

0.181 ±0.008 


The integral of in the measured range of 0.0025 < x < 0.7 is calculated using the data points. The 
contribution from the unmeasured region is extracted again from the fit. The various contributions are 
listed in Table [ 7 ] and the dependence of T/ s on the lower limit of the integral is shown in Fig. [ 8 ] The 
contribution of the measured x range to the integral corresponds to 93.8% of the full first moment, while 
the extrapolation to 0 and 1 amounts to 3.6% and 2.6%, respectively. Compared to the previous result |3j, 
the contribution of the extrapolation to x = 0 is now by about one third smaller than before due to the 
larger x range of the present data. The value of the integral for the full x range is 

rf s = 0.181 ±0.008 (stat.) ±0.014 (syst.) . (14) 




Figure 7: Values ofxg^ s (x) at Q 2 = 3 (GeV/c) 2 com- Figure 8: Values of f 1 as a function of x min . 

pared to the non-singlet NLO QCD fit using COM- The open circle at x = 0.7 is obtained from the fit. The 

PASS data only. The errors bars are statistical. The arrow on the left side shows the value for the full range, 

open square at lowest x is obtained with gj taken from 0 < x < 1. 
the NLO QCD fit. 

The uncertainty of Eft is dominated by the systematic uncertainties. The largest contribution stems 
from the uncertainty of the beam polarisation (5%); other contributions originate from uncertainties in 
the combined proton data, i.e. those of target polarisation, dilution factor and depolarisation factor. The 
uncertainties in the deuteron data have a smaller impact as the first moment of is smaller than that 
of the proton. The uncertainty due to the evolution to a common Q 2 is found to be negligible when 
varying Q 2 between 1 (GeV/c) 2 and 10 (GeV/ c ) 2 . The overall result agrees well with our earlier result 
= 0.190 ± 0.009 ±0.015 in g. 

The result for T^ s is used to evaluate the Bjorken sum rule with Eq. ( [To] ). Using the coefficient function 
Cf s (2 2 ) at NLO and a s = 0.337 at Q 2 = 3 (GeV/c) 2 , one obtains 
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|£a/#v| = 1.22±0.05 (stat.) ±0.10 (syst.). 


(15) 


The comparison of the value of \gh/gw\ from the present analysis and the one obtained from neutron /3 
decay, |#a/£v| = 1.2701 ±0.002 291, provides a validation of the Bjorken sum rule with an accuracy 
of 9%. Note that the contribution of Ay cancels in Eq. < [T2| ) and hence does not enter the Bjorken sum. 
Higher-order perturbative corrections are expected to increase slightly the result. By using the coefficient 
function C^ s at NNLO instead of NLO, |gA/gv| is found to be 1.25, closer to values stemming from the 
neutron weak decay. 


8 Conclusions 


The COMPASS Collaboration performed new measurements of the longitudinal double spin asymmetry 
Aj(jc,< 2 2 ) and the longitudinal spin structure function g\(x,Q 2 ) of the proton in the range 0.0025 < x < 
0.7 and in the DIS region, 1 < Q 2 < 190(GeV/c) 2 , thus extending the previously covered kinematic 
range [ 3 ] towards large values of Q 2 and small values of x. The new data improve the statistical precision 
of gj (x) by about a factor of two for x < 0.02. 


The world data for y' 1 , yj 1 and y“ were used to perform a NLO QCD analysis, including a detailed 
investigation of systematic effects. This analysis thus updates and supersedes the previous COMPASS 
QCD analysis [17], It was found that the contribution of quarks to the nucleon spin, A£, lies in the 
interval 0.26 and 0.36 at Q 1 = 3 (GeV/c) 2 , where the interval limits reflect mainly the large uncertainty 
in the determination of the gluon contribution. 


When combined with the previously published results on the deuteron [17], the new y^ data provide a 
new determination of the non-singlet spin structure function y^ s and a new evaluation of the Bjorken 
sum rule, which is validated to an accuracy of about 9 %. 
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Table A.l: Values of A ’j’ and gj as a function of x at the measured values of Q 2 . The first uncertainty is statistical, 
the second one systematic. 


x range 

(x) 

Q Z ((GeV/c) z ) 


S V \ 

0.0025 - 

0.004 

0.0035 

1.03 

0.059 ±0.029 ±0.014 

1.79±0.87±0.45 



0.0036 

1.10 

—0.004±0.027±0.012 

—0.12±0.81±0.37 



0.0038 

1.22 

0.002 ±0.032 ± 0.012 

0.05 ±0.98 ±0.37 

0.004 - 

0.005 

0.0044 

1.07 

0.006 ±0.021 ±0.008 

0.15±0.50±0.19 



0.0045 

1.24 

0.021 ±0.020±0.008 

0.53±0.51 ±0.20 



0.0046 

1.44 

0.023 ±0.022 ±0.011 

0.60 ±0.59 ±0.28 

0.005 - 

0.006 

0.0055 

1.11 

0.009 ±0.024 ±0.011 

0.18±0.46±0.21 



0.0055 

1.36 

0.026 ±0.020 ±0.008 

0.56±0.42±0.17 



0.0056 

1.68 

0.022 ±0.020 ±0.008 

0.51 ±0.47 ±0.18 

0.006 - 

0.008 

0.0069 

1.14 

0.033 ±0.020 ±0.009 

0.50±0.32±0.14 



0.0069 

1.50 

0.041 ±0.015 ±0.007 

0.71 ±0.27±0.12 



0.0071 

2.02 

0.006 ± 0.014 ± 0.007 

0.12±0.27±0.13 

0.008 - 

0.0010 

0.0089 

1.17 

0.007 ±0.027 ±0.013 

0.08±0.32±0.16 



0.0089 

1.62 

0.029 ±0.018 ±0.007 

0.40 ±0.25 ±0.10 



0.0090 

2.41 

0.015 ±0.014 ± 0.006 

0.24 ±0.23 ±0.09 

0.010 - 

0.014 

0.0116 

1.21 

0.044 ±0.026 ±0.013 

0.41 ±0.24±0.12 



0.0117 

1.75 

0.040 ±0.017 ±0.011 

0.42±0.18±0.11 



0.0120 

2.92 

0.044 ±0.011 ±0.005 

0.56±0.14±0.07 

0.014 - 

0.020 

0.0164 

1.26 

0.087 ±0.034±0.015 

0.58 ±0.22 ±0.11 



0.0165 

1.92 

0.100 ±0.020 ±0.011 

0.77±0.16±0.09 



0.0168 

3.74 

0.063 ±0.011 ±0.006 

0.60±0.10±0.06 

0.020 - 

0.030 

0.0239 

1.55 

0.072 ±0.030 ± 0.016 

0.36±0.15±0.08 



0.0240 

2.49 

0.079 ±0.025 ±0.011 

0.45±0.14±0.07 



0.0246 

5.16 

0.079 ±0.011 ±0.008 

0.545 ±0.077 ±0.061 

0.030 - 

0.040 

0.0341 

2.18 

0.103 ±0.035 ±0.016 

0.39±0.13±0.06 



0.0343 

3.50 

0.099 ±0.041 ±0.018 

0.43±0.18±0.08 



0.0347 

7.07 

0.083 ± 0.015 ±0.013 

0.421 ±0.075 ±0.066 

0.040 - 

0.060 

0.0473 

2.65 

0.128 ±0.040 ±0.023 

0.37±0.12±0.07 



0.0480 

5.00 

0.136 ±0.026 ±0.016 

0.449 ±0.086 ±0.057 



0.0492 

10.4 

0.103±0.015±0.011 

0.378 ±0.056 ±0.043 

0.060 - 

0.100 

0.0740 

4.91 

0.147 ±0.031 ± 0.016 

0.308 ±0.066 ±0.036 



0.0754 

10.7 

0.203 ±0.020 ± 0.017 

0.465 ±0.047 ±0.043 



0.0800 

19.7 

0.129 ±0.023 ± 0.021 

0.292 ±0.052 ±0.050 

0.100 - 

0.150 

0.119 

8.23 

0.291 ±0.038 ±0.024 

0.397 ±0.051 ±0.035 



0.121 

17.8 

0.263 ±0.028 ±0.021 

0.372 ±0.040 ±0.031 



0.125 

31.7 

0.242 ± 0.034 ± 0.021 

0.337 ±0.048 ±0.030 

0.150 - 

0.200 

0.171 

12.9 

0.299 ±0.045 ±0.027 

0.279 ±0.042 ±0.026 



0.172 

26.9 

0.316 ±0.045 ±0.036 

0.298 ±0.042 ±0.035 



0.175 

43.8 

0.344 ±0.050 ±0.029 

0.318 ±0.046 ±0.028 

0.200 - 

0.250 

0.222 

16.1 

0.405 ±0.060 ±0.043 

0.273 ±0.040 ±0.030 



0.222 

32.1 

0.340 ±0.066 ±0.035 

0.227 ±0.044 ±0.024 



0.224 

52.4 

0.268 ±0.060 ±0.045 

0.174 ±0.039 ±0.030 

0.250 - 

0.350 

0.289 

21.7 

0.397 ±0.057 ±0.035 

0.176 ±0.025 ±0.016 



0.290 

42.1 

0.374 ±0.077 ±0.050 

0.160 ±0.033 ±0.022 



0.296 

66.3 

0.392 ±0.062 ±0.035 

0.159 ±0.025 ±0.015 

0.350 - 

0.500 

0.403 

28.4 

0.396 ± 0.086 ± 0.051 

0.085 ±0.018 ±0.011 



0.405 

53.1 

0.40 ±0.12 ±0.06 

0.079 ±0.024 ±0.011 



0.413 

85.1 

0.631 ±0.088±0.054 

0.114 ±0.016 ±0.010 

0.500 - 

0.700 

0.561 

29.8 

0.42 ±0.23 ±0.10 

0.028 ±0.016 ±0.006 



0.567 

50.4 

0.75±0.23±0.12 

0.044 ±0.013 ±0.007 



0.575 

96.1 

0.87 ±0.15 ±0.09 

0.0429 ±0.0071 ±0.0048 
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